Zero temperature damping of Bose-Einstein condensate oscillations by 

vortex-antivortex pair creation 
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We investigate vortex-antivortex pair creation in a supersonically expanding and contracting 
quasi-2D Bose-Einstein condensate at zero temperature. For sufficiently large amplitude condensate 
oscillations, pair production provides the leading dissipation mechanism. The condensate oscilla- 
tions decay in a nonexponential fashion, and the dissipation rate depends strongly on the oscillation 
amplitude. These features allow to distinguish the decay due to pair creation from other possible 
damping mechanisms. Experimental observation of the predicted oscillation behavior of the super- 
fluid gas provides a direct confirmation of the hydrodynamical analogy of quantum electrodynamics 
and quantum vortex dynamics in two spatial dimensions. 
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The process of electron-positron pair creation is well- 
established in quantum electrodynamics since the sem- 
inal work of J. Schwinger pj. Later on, it became 
apparent that the hydrodynamics of vortices in two- 
dimensional (2D) superfluids can be mapped onto 2+ ID 
electrodynamics with vortices playing the role of charged 
particles, and phonons the role of photons [2j. In this 
analogy, the superfluid density and the supercurrent are 
acting as the magnetic and electric fields on the vortices 
whose circulation is the charge. The Schwinger vacuum 
breakdown is a phenomenon occuring whenever the elec- 
tric field exceeds the magnetic field (in cgs units), which 
corresponds in the analogy to the instability of a super- 
sonic flow with respect to the spontaneous creation of 
vortex-antivortex pairs from the superfluid vacuum. 

Vortices in Bose-Einstein condensates have been ob- 
served and studied experimentally intensely in the last 
couple of years, e.g., in 0, 0, la, Ifl ■ Here, we suggest 
an experiment in a quasi-2D Bose-condensed gas reveal- 
ing the existence of irreversible condensate dynamics at 
zero temperature as the result of the Schwinger pair cre- 
ation instability. To argue that vortex-antivortex pair 
creation is the dominant source of dissipation, we use the 
fact that a quasi-2D BEC in a time-dependent harmonic 
trap has a peculiar feature: There is a time dependent 
transformation (the so-called "scaling transformation"), 
using which the problem can be exactly solved, because 
the time dependence can effectively be removed from the 
Hamiltonian 0, 0, [^} . This scaling property also holds 
for linearized equations describing the evolution of small 
density and phase perturbations (Bogoliubov quasiparti- 
cle excitations) , propagating on top of the moving super- 
fluid. Therefore, initial perturbations cannot grow, and 
the instability mechanisms known from classical hydro- 
dynamics play no role 1 1 11] . This stability against pertur- 
bations implies that at very low temperatures, conden- 
sate oscillations are practically undamped (it has been 
measured that the quality factor Q > 2000 |l2|'). 

On the other hand, the superfluid velocity in the scal- 
ing solution grows linearly towards the condensate bor- 



der, while the local density decreases. Therefore the 
outer region of the cloud is always supersonic and, ac- 
cording to the Landau criterion for superfluidity, can host 
instabilities. Vortices are nonlinear excitations above the 
superfluid ground state, so that they are not protected by 
the scaling symmetry, which holds for linear excitations. 
Spontaneous vortex-antivortex pair creation, analogous 
to the Schwinger process, is an intrinsic instability mech- 
anism, and constitutes a source of dissipation already 
at zero temperature, without any need for a symmetry 
breaking external perturbation. 

In the following, we explicitly analyze the Schwinger 
instability of a supersonically expanding and contract- 
ing BEC in a time dependent quasi-2D harmonic trap. 
We show that for sufficiently large condensate oscillations 
vortex-antivortex pair production provides the dominant 
dissipation mechanism. Furthermore, the condensate os- 
cillations decay in a nonexponential fashion and the dis- 
sipation rate depends strongly on the oscillation ampli- 
tude. These features allow one to distinguish experimen- 
tally the decay due to pair creation from the previously 
studied damping mechanisms. We note that the sug- 
gested zero temperature damping mechanism is intrin- 
sically different from that discussed in |1(J |, where the 
dissipation is due to the energy transfer from the radial 
condensate motion to the longitudinal modes in an elon- 
gated cylindrically-symmetric condensate. This mecha- 
nism can only work if the condensate is sufficiently long, 
whereas we confine ourselves to the case of a quasi-2D 
sample, for which any motion along the z-axis is sup- 
pressed. 

The analogy of 2D vortex dynamics with electrody- 
namics is most easily established by noting that the ex- 
pression for the 2D Magnus force Fm = 27rpe z x (X — v s ) 
leads to the identification of E = pv s x e z and B = —pe z 
with the "electric" and "magnetic" fields, by comparing 
with the Lorentz force Fl = q(E + X x B). Here, X and 
v s are vortex and local superflow velocities, respectively, 
and p is the local density The circulation (27r in our units 
with h = m = 1) is the "charge" q (cf., e.g., [13, E3l ) 
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The self-energy of a (widely separated) single vortex pair 
is 2E® = 2npA, with A = \n(R/a c ), where R is the size 
of the pair. We will use in this expression for the pair 
energy that the vortex core size in a dilute superfluid is 
given by a c — l/c s , where c s is the speed of sound. The 
inertial rest mass of a vortex stemming from compress- 
ibility, m v — E®/c 2 , is (for large condensates) of "elec- 
trodynamical" origin: It stems from the self-interaction 
of a moving vortex with the long-range flow and den- 
sity fields it induces inside the surrounding superfluid 
medium. Since the "electromagnetic" fields (the density 
and velocity perturbations) represent "relativistic" par- 
ticles (phonons), the vortex mass diverges if the velocity 
of the accelerated vortex approaches the speed of sound, 
in the same manner in which the mass of a charged ul- 
trarelativistic particle diverges in conventional electrody- 
namics. We assume in what follows that other possible 
contributions to the vortex mass (see, e.g., the backflow 
mass contribution discussed in ^(j) remain regular if the 
local superfluid velocity approaches the speed of sound. 
These contributions are therefore subdominant for "rel- 
ativistically" moving vortices. 




FIG. 1: An oscillating, cylindrically symmetric quasi-2D con- 
densate. The shaded region designates the region of space in 
which the speed of sound is exceeded by the oscillating con- 
densate, and vortex pair creation takes place; H — H(t) is 
the horizon location and R = R(t) the Thomas-Fermi radius 
of the condensate. 

We consider a quasi-2D superfluid Bose gas in a 
time-dependent isotropic harmonic trapping potential 
V(x,t) — ^uj 2 (t)(x 2 +y 2 ), with x = (x,y). It is a well- 
established fact that the hydrodynamic solution for den- 
sity and velocity of motion in a harmonic potential with 
arbitrary time dependence may be obtained from a given 
initial solution by a scaling procedure 0, EJ • Defining the 
scaled coordinate vector r& = x/b, the rescaled density 
and velocity are given by 



P (x,t) = pa{n) = §r 1 1 



v s (x,t) = -X. 
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(1) 

(2) 



Here, we assume the superfluid to be described initially 
within the Thomas-Fermi (TF) approximation (that is, 
the condensate is large enough to neglect the quantum 
pressure); po is the initial central density and Rq the 



initial TF radius, such that R = R(t) = b(t)RQ is the 
instantaneous TF radius of the cloud. The energy func- 
tional has in the TF approximation the form 



d 2 n 




(3) 



with g the interaction strength, which depends in the 
present quasi-2D case on the tight confinement in z- 
direction and on the density of the condensate 0] . This 
leads to an effective Hamiltonian for the dynamical vari- 
able b, 



a l2 , a 2/A>2 , P 



£(b,b)=[-b 2 + -^(t)b 2 + ^ 



(4) 



where a — ttpoRq/6 and j3 — np^gRg/S. 

Consider a situation in which the external trap fre- 
quency is changed from oj- ln to w/ < Wj n , on a time 
scale much less than the inverse initial trap frequency. 
As a consequence, the gas undergoes a large amplitude 
monopole oscillation with frequency 2w/ yjj. At suffi- 
ciently low temperatures (below the Kosterlitz-Thouless 
temperature), the initial state of the superfluid con- 
tains bound vortex-antivortex pairs, i.e. topological exci- 
tations, which can be unbound by the action of the Mag- 
nus force in the (time dependent) supersonic flow region. 
Indeed, for an oscillating condensate, there exists a re- 
gion, the border of which is called horizon (cf. Fig.QJ, 
where the superfluid velocity magnitude v s is larger than 
the local sound velocity c s — ^fgp. The speed of sound 
is exceeded at the horizon radius 



H(t) = 



b(t)R 
^7 2 W + 1 



(5) 



where 7 = \/2bb/ u>m- 

Beyond the horizon, the vortices and antivortices get 
accelerated during condensate evolution and separate at 
local superflow velocities larger than that of sound. This 
is analogous to the Schwinger pair creation process in 
quantum electrodynamics. It is important to recognize 
that the flow we consider is inhomogeneous and time de- 
pendent by default. Consequently, the argument that 
there is no pair creation possible because one could al- 
ways use the underlying Galilean invariance to "trans- 
form away" the background flow, does not apply to our 
situation. 

In a simple model of the 2+1D vacuum pair creation in- 
stability, which exploits directly the analogy to Schwinger 
pair creation in quantum electrodynamics, the pai r pro- 
duction rate T per unit area can be written as |19| 
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where we have defined T — E 2 c 2 — B 2 c\ and set, within 
logarithmic accuracy, the vortex pair size in E® equal 
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to the Thomas-Fermi radius of the condensate. The 
above relation holds for locally supersonic motion, i.e., 
ii \E\/\B\ > c s {T > 0). The value of the prefactor in 
front of the exponential in the above expression is sub- 
ject to changes which are due to the microscopic details 
of vortex motion. We display its value, stemming from 
taking literal the analogy to quantum electrodynamics 
also on the level of quantum fluctuations (to one loop or- 
der), for numerical concreteness. However, the behavior 
of r for > c s is dominated by the hydrodynami- 

cal exponent, whose value is independent of microscopic 
physics, and specifically by the n — 1 term in the above 
sum. 

Assuming that the vortex density is low, the energy 
dissipation rate, e, is obtained by multiplying Eq. © by 
the rest energy of the widely separated vortices 2E®, and 
integrating over the area of the TF domain. This results 
in 
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( 7 2 + 1 )13/4 I g 



(7) 



where we introduced the function 



F(X) 



E- 



,3/2 



(8) 

Since 7 is proportional to bb, the Schwinger dissipa- 
tion rate @ can give rise to a measurable effect only if 
the condensate oscillation amplitude is sufficiently large, 
which implies u>f -C Lo- m ■ In order to provide some analyt- 
ical results, we consider a simple quasistationary pertur- 
bation theory approach. Consequently, we assume that 
the dissipation rate is small and therefore that the en- 
ergy of the system in Eq. J3J is a slowly varying function 
within each oscillation period. Then the equation of mo- 
tion for the scaling parameter b can be found from 



dt 



£(b, b) 



(9) 



In the absence of dissipation (e = 0), the range of b is 
between b mnl — 1 and 6 max = u)- m /u)f. Since b Illm = 1 <C 
frmax, we can approximately set 6 m i n ~ 0. One can then 
write 7 2 = 2(w 2 /w 1 2 n )6 2 (& 2 nax — b 2 ). In a dilute gas, in the 

TF limit, the argument of F is large, A 2 ^S> g\/"f 2 + 1, 
and the dynamical equation (0 for b takes the simpler 
form: 



b + uj 2 f b - 



V 



7 hi 



b 3 



-b'b 



(10) 



lost in a period is then given by 



Ik = 



7tt 

1024" 



V^b 



f Ll2 



(12) 



The energy decrease rate for b n 
equation 



is obtained from the 



-j;o{o nmx ) — — 1e- 
at 7r 



(13) 



Thus one obtains for the oscillation peak value the fol- 
lowing expression 



c(t) 



c(0) 



(1 



Vb 10 

"max 



(0)w/t) 



1/10 ' 



(14) 



where V — -^(u>f /u>i n ) T>. Our perturbation theory 
approach is valid as long as b^ ax (0)V <C 1. 

The damping of condensate oscillations due to vortex- 
antivortex pair creation is represented in Fig.|2 where we 
show the numerical solution of the dynamical equation 
Eq. IjlOl) (grey solid line), the approximate solution for 
the peak amplitude (black solid line) , and for comparison 
the free oscillation without pair creation (dashed line). 
The parameters used in the numerical integration for the 
plot are N — 10 4 , g = 1, and for the final trapping 
frequency LUf — 0.1 u>i n . These parameters are consistent 
with the argument of F being large, A 2 3> g-y/7 2 + 1, so 
that Eqs. i|TU|) - ((T^J) hold. The envelope 6 max (i) is seen 
to decay very slowly and in a nonexponential fashion, 
governed by the TF exponent in Eq. 1)14(1. For realistic 




where the constant 



V 



(in [lJg~N~fc 



(-l) nH 



(11) 



Using the equation (|10|l , bb can be expressed in terms of 
b only, 6 2 6 2 = uj 2 b 2 (b 2 lax — b 2 ). The oscillation energy 



FIG. 2: Damping of condensate oscillations due to vortex- 
antivortex pair creation, with N = 10 4 , g = 1, and u)j = 
0.1 t^in, where the radius 7? is in units of the original Thomas- 
Fermi size Ro- The black solid line is the envelope b max from 
Eq. ill II . The grey solid line is the damped breathing mode 
oscillation obtained from numerically solving Eq. 1101 . For 
comparison, the dashed line represents the oscillation of the 
superfluid gas without pair creation taking place. 
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parameters, we conclude from Fig. [5] that an observable 
damping effect for the condensate oscillations is obtained. 

The scaling parameter evolution can be described by 
Eq. J5J only for sufficiently short times, when the total 
density of vortices produced is still low. At later times, 
the vortex-antivortex plasma can decrease the superfluid 
current in the same way as the electron-positron plasma 
can screen the electric field. This is an interesting col- 
lective effect, which however requires a more elaborate 
treatment. 

We described an intrinsic damping mechanism for large 
amplitude condensate oscillations in a quasi-2D Bose gas 
at zero temperature. The dissipation originates from 
spontaneous creation of vortex-antivortex pairs and de- 
pends on the shape and dynamics of the supersonic 
flow region. The results we presented therefore depend 
strongly on the oscillation amplitude. This feature can 
be used to distinguish the effects of pair production from 
other possible dissipation mechanisms. The scaling solu- 
tion not only exists for the discussed monopole modes, 
but also for quadrupole oscillations, so that, e.g., effects 
resulting from a rotating superfluid on the pair creation 
process may be studied. Observation of the predicted 
oscillation behavior of the superfluid gas provides direct 



confirmation of the hydrodynamical analogy of quantum 
electrodynamics and quantum vortex dynamics in two 
spatial dimensions, and would put this analogy to its 
first real experimental test. Such confirmation would, 
then, give further motivation to the program of studying 
analogies between high energy physics, cosmology and 
condensed matter systems [2£|. 

The outlined mechanism for dissipation is not confined 
to quasi-2D samples. In strongly elongated 3D conden- 
sates, the scaling solution also applies, and the vorticity is 
generated in the form of vortex rings, with the total vor- 
ticity integrated over the sample volume still zero. How- 
ever, as already mentioned above, for a 3D condensate 
the effect of vortex ring creation can be masked by pos- 
sibly stronger damping mechanisms, like the parametric 
resonance discussed in [Tof . 
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